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The variational principle, differential equations, and boundary 
conditions governing the cross-sectional distortions due to inertia 
loading of a two-dimensional model of a thin monocoque wing are shown. 
A theoretical analysis of this simplified model is made in order to 
determine the nature of the coupling between the cross-sectional modes 
and the spanwise deformation modes. General solutions are obtained in 
finite-difference form for arbitrary cross sections and an exact solu- 
tion is presented for a parabolic-arc cross section of constant cover 
thickness. The application of these results in evaluating the coupled 
frequencies of the actual structure is discussed. Frequencies evalu- 
ated for a parabolic-arc monocoque beam show good agreement with 
experimental values. 


INTRODUCTION 


Chordwise modes - that is, modes whose ehordwise variation in 
deflection is such that more than one node line exists at every cross 
section along the span - have been observed in several dynamic inves- 
tigations of pure monocoque and multiweb wing models. For example, in 
cm extensive series of tests of aero dynamically heated multiweb wing 
structures the principal mode of failure was established to be chord- 
wise flutter. (See ref. 1.) Also, in an investigation of the vibra- 
tional characteristics of circular-arc monocoque beams, discussed in 
reference 2, several "extra” resonant conditions were observed experi- 
mentally which, although referred to as panel modes, are more accu- 
rately described as chordwise or cross-sectional modes. In addition, 
the results of reference 2 indicate that the accuracy of the predicted 
overall bending and torsion frequencies, even after the inclusion of 
important secondaxy effects, such as transverse shear and shear lag 
in bending and restraint of warping in torsion, are considerably 
impaired by the presence of these chordwise modes. For instance, the 



2 


* 


calculated frequency of even the first symmetrical torsion mode (mode 2 
in table II of ref. 2) differed from the one obtained experimentally by 
72 percent while the third symmetrical mode differed by more than 
250 percent. 

The influence of cross-sectional distortions on the overall beam 
modes was investigated in references 3 and ^ • These investigations, 
however, were not concerned primarily with the prediction of cross- 
sectional modes but rather with the influence of cross-sectional flexi- 
bilities on the overall beam modes. Furthemore, neither analysis is 
applicable to monocoque or multiweb wings with thin airfoil-type cross 
sections such as those used in references 1 and 2. 

The purpose of the present paper is to investigate the chordwise 
modes and frequencies of beams with thin airfoil-type cross sections 
and to determine the influence of cross-sectional flexibilities on the 
overall beam frequencies. A simplified model which embodies the prin- 
cipal mechanisms of cross-sectional and sparwise flexibility is analyzed. 
The differential equations and boundary conditions for the simplified 
model appropriate for wings of arbitrary airfoil-type cross section are 
derived. An exact solution is found for the model for the particular 
case of a monocoque wing of parabolic cross section. In addition, a 
general solution is presented in finite-difference form by applying the 
calculus of variations directly to the energy expressions. A discussion 
is made of the use of these results to predict the frequencies of both 
the coupled cross-sectional modes and the overall beam modes of a mono- 
coque wing presented in reference 2. 


SYMBOLS 


c chord length of beam 

E modulus of elasticity 

h distance from x-axis to median line of skin 

ho maxi mum value of h 

I moment of inertia about median lire of skin, t^/l2 

k-^k^k spring constants 

M 


N 


moment (see fig. 2) 

number of discrete stations in finite-difference analysis 


h! H 4- 



radius of curvature of midplane of beam 

tangential coordinate 

maximum kinetic energy 

skin thickness 

maximum strain energy 

displacement of midplane in x-direction 
displacement in s-direction, h’w + u 
shear (see fig. 2) 

deflection in z-direction (positive upward) 
deflection normal to s-direction 
chordwise coordinate 
spanwise coordinate 

coordinate perpendicular to x and y coordinates 
frequency coefficient 

frequency coefficient for uncoupled chordwise mode 
strain in s-direction 
strain in x-direction 

nondimen sional distance "between stations on the half- chord 
normal midplane force in the beam 

mass density of cover 
natural frequency 

uncoupled cross-sectional frequency 



k 


®e,l 

u>exper 

cu r, B 
o>r,T 


fundamental or first uncoupled cross-sectional frequency 
experimental natural frequency 
rigid-body bending frequency, 




k 

cpt 


^k 

rigid-body torsion frequency, \ — — 

cpt 


Primes on symbols denote differentiation with respect to x; primes 
on matrices denote the transpose of a matrix. 
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Idealization 

A rigorous approach to the problem of the vibrational characteris- 
tics of a monocoque wing which includes both spanwise and cross-sectional 
flexibilities . would, of necessity, involve some type of shell analysis. 

A more appealing although less rigorous approach, however, is to investi- 
gate first the uncoupled chordwise and spanwise characteristics of the 
structure independently and then determine the influence on the natural 
frequencies of the coupling of these motions. This can be done in a 
manner similar to that used in references 3 and 4 by mounting a typical 
cross section of the wing on elastic supports in such a way as to allow 
for the necessary overall degrees of freedom. 

The only structures considered in this paper have uniform properties 
in the spanwise direction with a symmetrical curved airfoil cross section, 
as shown in figure l(a). The wing or beam is of either pure monocoque or 
multiweb construction with no, or relatively few, bulkheads. An idealized 
structure that can be easily analyzed but still possesses essentially the 
same cross-sectional and spanwise flexibilities as the original is shown 
in figure l(b). 

The idealized structure consists of two curved beams attached at 
the ends to elastic supports in such a manner as to maintain the angle 
between the beams at each end. The idealization is the same for both 
the monocoque and multiweb prototypes. Consequently, the spanwise webs 
of the multiweb structure are assumed to contribute nothing to the chord- 
wise stiffness of the wing. The webs, however, do tend to prevent any 
relative motion between the top and bottom covers. Hence the additional 
assumption is made that deflections of the top and bottom beams in the 
idealized structure are the same. The depth of the individual beams in 
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the idealized structure is taken to be the same as the thickness of the 
covers in the prototype. The stiffness of the elastic supports or springs 
is chosen to duplicate the bending or torsional stiffness of the proto- 
type, depending on whether the motion is symmetrical or anti symmetrical 
about the midchord. 
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Method of Analysis 

The structure to be analyzed is shown in figure 2 along with the 
coordinate system and sign convention to be used. As the structure is 
symmetrical about the x-axis and the distance between the top and bottom 
segments is maintained, only half of the structure need be considered. 
Each segment is assumed to carry midplane forces in addition to trans- 
verse shear and bending moments (see fig. 2), and will be considered to 
act in a manner consistent with elementary beam theory (i.e., plane sec- 
tions remain plane). 

The strain in the midplane of the beam is given by 


€ 


S " 


w 

r 


+ 


du 

ds 


(i) 


where 

v deflection normal to s-direction 

r radius of curvature of midplane 

u displacement in s-direction 

s tangential coordinate 

Since the cross section of the prototype being analyzed is that of 
a thin airfoil section and deformations are considered small, the tan- 
gential coordinate s can be approximated by the chordwise coordinate x, 
and the curvature by -h ff . Furthermore, the deflection w can be con- 
sidered to be the same as the deflection in the z-direction, and the dis- 
placement in the s-direction can be written as 

u = h'w + u (2) 

where w is the deflection in the z-direction and u is the displace- 
ment in the x- direction. Hence the expression for strain in equation (l) 
becomes 


e x = k' wt + u * 


( 5 ) 
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If this expression is used for the strain, and w and u are considered 
as amplitudes of displacement, the maximum strain energy for one of the 
covers of the idealized structure vibrating in a natural mode becomes 

U = i f El(v") 2 dx + f A(h'w' + u' )dx + ~ ^w^c) + i kpw 2 (0) (4) 

d J 0 ^0 ■ 4 

where 

El bending stiffness of beam L 

1 

A normal midplane force 4 

4 

k^,k 2 spring constants 4 

The first term in equation (4) represents the contribution to the strain 

energy due to bending of the idealized structure, the second represents 

the strain energy due to the midplane force, and the third and fourth 

terms are the potential energy of the spring mounts. • 

In the present analysis the force A will be considered constant 
and equal to the average midplane force; hence, with the use of equa- 
tion (3), 



+ u' )dx 


(5) 


Furthermore, if inextensional theory is assumed u' vanishes in 
equation (4), and equation ( 5 ) becomes 


A 


1 

c 



Eth'v'dx 


(6) 


The maximum kinetic energy of the idealised structure vibrating in 
a natural mode is 


T 




jWt^cLx 


(7) 


where o> is the natural circular frequency of the mode un der considera- 
tion and p is the mass density of the cover. 
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A natural mode of vibration must satisfy the variational equation 

5(U - T) = 0 (8) 

where the variation is taken independently with respect to w. Applica- 
tion of equation (8) and the usual methods of the calculus of variations 
results in the following differential equation and boundary conditions: 

(EIv" )" - Ah" - ptu) 2 w = 0 (9) 


(EIw" ) ' 


(EIw") ’ 


(EIw" )S w' 


= 0 


X=C 


- Ah’ - w 

2 


6w 


= 0 


x=c 


(EIw" ) Sw 1 


I x=0 


= 0 


k 2 

Ah + — w 


5w I 


x=0 


= 0 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


where A is given by equation (6). 

If the cover thicknesses of the original wing are constant and if 
the half-depth h of the wing is approximately parabolic (as, for 
example, the circular-arc airfoils used in refs. 1 and 2) an exact solu- 
tion to the differential equation is possible and is discussed in the 
next section. However, for an airfoil of arbitrary shape and variable 
cover thickness it is unlikely that an exact solution can be found. 

For these cases the use of some approximate variational procedure is 
advisable. One such approach is presented in the appendix. Here, by a 
procedure similar to that discussed in reference 5> the energy is 
expressed in matrix form in terms of the deflections at discrete points 
by means of finite- difference relationships. The variational procedure 
results in a frequency equation in matrix form from which the mode 
shapes and frequencies of the structure can be calculated. 


Exact Solution of Idealization of the Monocoque Be.am 
of Parabolic- Arc Cross Section 

For an airfoil of doubly symmetrical cross section and constant 
cover thickness whose half- depth h can be represented by a parabolic 
equation, an exact solution to the differential equation for the cross 
section mounted on springs is possible. In order to take advantage of 
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the double symmetry the origin of the coordinate system is now relocated 
at the midchord of the cross section and the nondimensional coordinate 
2x 

| = — — 1 is used* Thus, the half-depth can be represented as 


h = h Q (l - I 2 ) 

Furthermore, the spring constants are identical, so that k-^ = kQ 
Equation (9) can now be written in nondimensional form as 

d^w Ac 2 d 2 h 4 n 

d| 4 4EI dg 2 


(14) 

= k. 


(15) 


where 


4+2 
4 captor 

a = 

16EI 


(16) 
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A = 


4h 0 tE 


w(-l) 



w(|)d| 


(IT) 


Because of the constant cover thicknesses of the original wing, the 
more stringent conditions of inextensional theory are not necessary in 
determining X. It need only be assumed that there is no overall short- 
ening or stretching of the cross section; then the integral of u 1 will 
vanish in equations (4) and (5) instead of u* being zero everywhere. 


Because, in the general solution to the differential equation, the 
symmetrical and antisymmetrical parts are easily recognizable, they are 
presented separately. 


Symmetrical solution *- For the symmetrical solution the boundary 
conditions, as given by equations (lO) and (ll), can be written in non- 
dimensional form as 


d| 2 


1=1 


0 


(18) 


and 



d^ w 

d| 


4 V , k 

5 + — t - 2 V 


= 0 


1=1 


(19) 
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The symmetrical solution to the homogeneous part of equation ( 15 ) is 

w = Aq cos <x£ + Bq cosh o| (20) 

and the particular solution is simply 


w = 


Ah Q c 
2EI0 4 


( 21 ) 


Thus 


w = Aq cos a{ + Bq cosh a| + 


AhoC^ 

2EIO* 


( 22 ) 


Substitution of equation (22) into the boundary conditions for the sym- 
metrical solution (eqs. (l8) and ( 19 )) yields, respectively, 


a _ cosh a ^ 
Aq Bq 

cos a 


(23) 


B i 

— cosh a - 
A 


^ h 0fq _ _3kcl_) 
C \ 4EtW 


El ^2L_(tan a + tanh a) - k 

C 5 


(24) 


Note that three unknowns appear in equations ( 23 ) and (24): Aq, 

Bq, and A. However, by substituting equation (22) for w into equa- 
tion (17) for A and performing the necessary integration, a third 
expression is obtained: 


8h tE / 

A = — ^ — Bq cosh a^-2 


tan a + tanh a' 


(25) 


If Bi 


cosh a 


in equation ( 25 ) is replaced by its identity in equa- 

A 

tion (24), the characteristic or frequency equation is found, and can 
be written as 


1 = 


W (t) [a* - <t« a + 1»oh aj) [l - 

3 

a*(tan a + tanh a) - 


( 26 ) 
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Substitution of equations ( 23 ) and ( 25 ) into equation (22) yields the 
mode shape 


cos cosh a| 

w = + 


cos a 


cosh a 


tan a + 

tanh a 

- 5 \( c ? 

a 


2EcMw 


3^ 

/c\3 

-L 

4Ea J+ 

\t) 


(27) 


Anti symmetrical solution .- In the anti symmetrical solution the mid- 
plane force A is zero, as can be seen from equation ( 6 ). In equation ( 6 ) 
the integrand of the integral expression for K is an anti symmetrical 
function about c/ 2 ; consequently its integral from 0 to c is zero. 

The differential equation for the cross section, given by equation (l5)> 
and the boundary conditions given by equations (lO) and (ll) therefore 
reduce to 


- o>v = 0 

dl 4 


( 28 ) 
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and 


a g w 

di 2 


= 0 


(29) 


l=i 


EI — - - v 

c ' d |5 2 


1=1 


= 0 


(30) 


Equations (28) to ( 50 ) are the differential equation and boundary 
conditions for a beam supported by springs at £ = -1 and | = 1. The 
solution is presented in most vibration text books, but for the sake of 
completeness it will be repeated here. For the antisymmetrical solution 


w 


A 2 sin a£ + B 2 sinh a£ 

Substitution of equation ( 31 ) into the boundary conditions yields 


(31) 


a _ sinh a ^ 
A p bp 

^ sin a ^ 


( 32 ) 


El(-A2 cos a + B2 cosh a) = ^^2 sin a + B2 sinh a) 


(33) 
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Thus, the characteristic or frequency equation becomes 


coth a - cot a 
2a 


3k 

4Ea^ 


(!) : 


and the mode shape is 


sin a4 j sinh a£ 
sin a sinh a 


(3M 


(35) 


NUMERICAL RESULTS FOR IDEALIZATION OF MONOCOQUE WING 
WITH SYMMETRICAL PARABOLIC CROSS SECTION 


Symmetrical Modes 


The exact frequency and mode- shape equations of the idealized struc- 
ture were obtained in the previous section. In equation (26) the fre- 


quency parameter a I where a> = a' 


- „2.^E \ 


nw 


of the symmetrical modes was 


The motion 


given as a function of t/h Q and the parameter 

of the idealized structure can be considered to be the result of the 
coupling of two types of modes of vibration: an infinite set of 

uncoupled cross-sectional modes which are obtained by omitting the 
springs, and a rigid-body mode represented by a rigid cross section 
mounted on springs. Consequently, if the vibrations are considered in 

Kf ) 3 


this light a more useful parameter than 


would be ux 


•,b/' 


'CD. 


where 


“>r,B 


by 


is the rigid-body translation frequency and is defined 

—7—, so that 
ptc 



Thus the frequency equation (26) becomes 



/ , v p k 
(±\ _ 1 

2a 


VV 48 j 

tan a + tanh a 


", . . 2 4 1 



(±\ £L + 1 

2a 


UJ 48 J 

tan a + tanh a 



(36) 


( 37 ) 
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and the corresponding mode shape (eq. ( 27 )) is given by 


tan a + tanli 


cos ai cosh a£ 
w = — f 

cos a cosh a 




°V,B 


The coupled frequency coefficient a was calculated from equa- 

/ t \ 2 /“ r b \ 2 

tion (37) for various values of f — J and ^ J and the results 

are shown in figure 3 . Here the frequency coefficient a for the first 

2 / t \ 2 

three modes is plotted as a function of ( - - --- \ for a range of (r— ) 


from 0 to 0 . 09 • As 


becomes greater than 1, all curves approach 


horizontal asymptotes whose values depend on — • These asymptotes have 

" ho 

the value of the uncoupled cross-sectional frequencies obtained by 
letting k = 0 which correspond to the case where no coupling exists 
between rigid and cross-sectional motion. Consequently the mode shapes 
associated with these frequencies are ones which are predominantly cross 
sectional in character - that is, they contain two or more nodal points. 

For k equal to zero, equation ( 37 ) reduces to the frequency equa- 
tion for the uncoupled cross-sectional mode: 


— 1 

4\tan a c + tanh ol c 


and equation ( 38 ) reduces to the corresponding uncoupled cross-sectional 
mode shape : 

cos a c | cosh a c | 2 , )jM 


cosh a„ 


Qc 1 * /t \ 2 + 

48 \hj 


where 



L 

1 

4 

4 

4 


( 41 ) 
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In figures 4 (a) and 4 (b), a c is shown plotted as a function of 



for the first two symmetrical frequencies. 


These plots serve to 


determine the uncoupled chordwise frequency which corresponds to the 
pure chordwise modes of an infinitely long wing. 


Values at 


O) 


01 


r, B 


=0 in figure 3 describe the condition where 


the ends of the cross section are pinned; that is, k = Hence, mode 


shapes for 






*r, B 


« 1 are also of predominantly chordwise character. 


For values of 


CD 


OD. 


V,B 


approximately equal to 1 the curves approach a 


vertical asymptote = 1 and correspond to the condition where 

the mode shapes are of the rigid-body type; that is, they contain no 

nodal points. As can be seen from figure 3 j when ( — 0 the curves 

\ h o/ 

fall on the asymptotes and no coupling exists. This uncoupling of the 
two motions could also have been surmised by examining the mode shape 

as given by equation (38). Note that for (— j = 0 the tip deflec- 
tions of the cross section in equation (38) are zero, and hence the 
stiffness of the springs can have no influence on the overall motion. 
The differences between the curves and their corresponding asymptotes 


for the other values of 


-V 


represent the change in frequency due 


to coupling between the chordwise and rigid motions. As can be seen 
from the plots, the greatest effect is felt when co < 


One drawback with the curves presented in figure 3 is the difficulty 
encountered in using these curves to calculate the coupled frequency, 
since the unknown coupled frequency coefficient a, appears in both the 
abscissa and ordinate of the curves. Thus the results of figure 3 have 
been replotted in figures 5(a) and 5(b) in terms of slightly different 
parameters. 


In figures 5 (a) and 5 (l) the ratio of the coupled frequency ad to 
the rigid-body frequency 0) r ^ is plotted as a function of the ratio 

of the fundamental uncoupled chordwise frequency cjd to the rigid- 

c , 1 

a ) , 

body frequency o> r fi . Figure 5 (a) presents the data for values of — — 

r, B 
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from 0 to li .8 and from 0 to 3-2. Figure 5(b) is simply an 

-n 

>B “cl 

enlargement of figure 5 (a) in the region where goes from 0 to 2 . 0 . 

r, B 

The values of to in the abscissa are based on the first uncoupled 

C , JL 

symmetrical mode for k = 0. Hence the abscissa depends only on the 
stiffness properties of the structure and not on the frequency parameter. 
The vertical asymptote of figure 3 appears in figure 5 as a horizontal 

line at “■■■ = 1 . 0 , while the horizontal asymptotes of figure 3 tranS- 
r, B 

form into straight lines which emanate from rhe origin. Note that only 
2 

the (— ) = 0 asymptotes appear in figure 5, since the asymptotes for 

Vho / 9 

2 /t v 2 

all other values of lie in such close proximity to the = 0 

asymptotes that they cannot be easily distinguished. 

Figure 5 shows that for values of q greater than g the 

influence of coupling is appreciable only for the first or lowest mode, 
on 

while for — - — < 1 all inodes investigated show substantial effects, 
r, B 

, Ph 0 frequencies of the coupled modes are easily obtained from figure 5^ 
since the abscissa depends only on the physical properties which, of 
course, are assumed to be known. The difficulty lies in establishing 
the characteristics of the mode shapes - that is, in determining whether 
the resulting mode is predominantly a rigid mode with no nodal points or 
predominantly a chordwise mode with two or more nodal points. 


In order to try to resolve this problem, the first three modes for 


( 


2 -I 

_t_\ = 0.018 and = 0.987 and 0.409 were calculated and are 

hj “r,B 

shown in figures 6 (a) and 6 (b). From figure 6 (a) it can be seen that 
the first mode is predominantly rigid-body irotion while the remaining 
modes consist of chordwise-type motions. Ir. figure 6 (b), however, the 

second mode is the rigid-body- type mode. Presumably, as 




®r,B 


is 


reduced still further the third mode becomes the rigid-body mode and so 
on. In order to illustrate the change in mcde shape due to coupling of 
chordwise and rigid-body motions, a sketch cf the second mode for 
,2 

( —\ = 0 . 0 l 8 and various values of 

UJ 

is shown in figure 6 (c). 


O) l 

— - — ranging from 0.026 to 0 . 9&7 

(» -r. 


r, B 


L 

1 
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4 

4 



H -=t -=t ^ 
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From the foregoing discussion it 
rigid mode is always the lowest mode. 


CD 


I -1 

can be seen that for — 1 > 1 the 


a ). 


o> 


For 


c,l 


r,B 


< 1 it is advisable to 


1 y n 

calculate the mode shape in order to determine which mode is actually 
the rigid-body mode. 


In reference k the behavior of rectangular bents with flexible 
members, mounted on either deflectional or rotational springs, has been 
analyzed In detail. 


It Is of interest to note that the assumption of h^ = 0 



in the present paper corresponds to the limiting case in reference k 
where the depth of the web is assumed to become zero while the thickness 
becomes infinite. For such a configuration the frequency equations of 
the two papers are found to be identical. This Is true for both the 
symmetrical and the antisymmetrical modes. 


Antisymmetrical Modes 


For the antisymmetrical modes, the frequency parameter a as given 

in equation (34) is a function of only. This lack of dependence 

of frequency on camber in the antisymmetrical modes is a consequence of 
the vanishing of the midplane force A. Therefore, the equations for 
the mode shapes and frequencies are those of the antisymmetrical vibra- 
tions of a beam mounted on spring supports. As was true for the symmet- 
rical modes, the motion of the idealized structure can be considered as 
a result of the coupling between cross-sectional modes and rigid-body 
modes. In the antisymmetrical case, however, the rigid-body mode is a 


pitching or torsional mode. 

rta 
is 




in which T 


Hence, a more useful parameter than 
is now defined by 



3 



3k 

ptc 


( 42 ) 


The parameter ca T can be looked upon as a rigid-body torsion mode 
represented by the well-known equation 
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where 



and 

T - Ptc5 

Ip = T" 


( 43 ) 

(44) 

(45) 


Note that Ip is the approximate polar moment of inertia for a thin 

cross section. The characteristic or frequency equation (eq. (34)) 
then be written as 

2 

3 coth a - cot a _ f 0 * r, t \ 

2 a \ a) / 


can 

(46) 


The mode shape is not an explicit function of the spring constant and 
remains 


sin at sinh ag 

sin a sinh a 


(47) 


The coupled frequency coefficient a was calculated from equa- 

2 

tion (46) for various values of f ] and the results are shown in 

“r,T/ 


1 ' 

figure 7. In figure 7 the first three modes of a are plotted as a 
, ... ^2 

function of 




The uncoupled cross-sectional modes are repre- 


sented by dashed horizontal lines, at a = 3-927 and 7- 0^9* ** or ^ = 0. 
The pin-ended condition for k = oo is represented by the dashed lines 
at a = n, 2*, and 3*. The rigid-body torsion mode is the dashed verti- 

O) 


cal line at 


= 1. 


The coupled modes are shown by the curved lines. 


(JD 


o> 


as a function of 


c,l 


is shown in 


®r,T 

a) in the abscissa are based on the first 

c, 1 


®r,T 

The corresponding plot of 

figure 8. The values of 
uncoupled antisyrametrical mode for k = 0. Hie results are plotted in 
a manner similar to that for the symmetrical modes but only one curve 
is presented for each mode, the results being independent of the 

to 


L 

1 

4 
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or camber effect. The dashed horizontal line at 


=1.0 represents 



the rigid-body torsion mode, and the dashed straight lines inclined to 
the abscissa represent the uncoupled chordwise modes. 

Contrary to the case of the symmetrical modes, the predominantly 
rigid-body mode with no nodal points is easily distinguished, as it is 
always the first or fundamental mode. 


APPLICATION OF METHOD 


The present paper has two objectives: (l) to analyze the chordwise 

vibrational characteristics of thin wings with airfoil-type cross sec- 
tions and (2) to obtain the effect of the coupling between the chordwise 
and spanwise motions of the structure. Both of these objectives can be 
accomplished by utilizing the results obtained from the analyses of the 
simplified model. 

In principle, the method presented is similar to that discussed in 
references 3 and 4. In the present investigation, however, emphasis is 
placed on the determination of the predominantly chordwise modes and 
frequencies of the structure in addition to the modes and frequencies 
of the overall beam motions. The results obtained from the analysis for 
the model with springs can be used to give both the chordwise and the 
spanwise coupled frequencies provided, of course, that the appropriate 
choice of the spring constant is made. The uncoupled cross-sectional 
modes and frequencies can be obtained from the analysis by letting the 
spring constant vanish. 

Since the rigid-body motion of the idealized structure depends only 
on the spring constant while the relative chordwise motion depends on 
the flexibility of its members, the value of the spring constant should 
be chosen so that the uncoupled rigid-body frequency (either transla- 
tional or torsional) of the idealized structure is the same as the span- 
wise, frequency of the prototype. For the symmetrical chordwise modes 
the spring constant should be determined from the frequencies of the 
spanwise bending modes by means of the equation 

k = ^ptc (48) 

while for the antisymmetrical chordwise modes the spring constant should 
be determined from the spanwise torsional frequencies by means of the 
equation 
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The spanwise frequency of the prototype should be obtained by the most 
appropriate method available, with the provision that the method be based 
on the assumption that the cross section of 'the structure remain 
undistorted. 

At this time a note of caution as to the application of the method 
presented is in order. The use of the idealized structure embodies the 
assumption that all cross sections distort identically. Hence the influ- 
ence of any bulkhead, rib, or other chordwise restraint is not included. 
Consequently, the method is applicable only to wings in which the bulk- 
head spacing is not very small. L 

1 

In order to demonstrate the use of the method, both the coupled 4 

chordwise and coupled spanwise modes and frequencies of a circular-arc 4 

monocoque beam used in the experimental investigation of reference 2 will 4 

be calculated. For the sake of completeness a short description of the 
specimen used and the results obtained will be presented in this paper. 

The aluminum monocoque beam was 140 inches long and 40 inches wide. 

The nominal dimensions and average radius of curvature are shown in 

h 

figure 9* The beam had a ratio — of approximately 0.04 and a ratio 

/ 1 \ 2 c 

J of 0.06735* The cross sections were maintained by steel bulk- 
heads located 10 inches from each spanwise tip. Eighteen bending and 
torsional modes were found, and in addition a number of so-called "panel" 
modes were located. (See figs. 5 and 7 of ref. 2.) 

From a comparison of the spanwise node lines shown in reference 2 
and the mode shapes shown in figure 6, it can be concluded that many of 
these panel modes can be described as chordwise modes. The variation 
in position of the chordwise node lines is the result of coupling of 
the symmetrical chordwise modes with the beam bending modes or coupling 
of the anti symmetrical chordwise modes with the beam torsional modes. 

Consequently, a family of chordwise modes is associated with each of 
the spanwise modes. An additional node line appears at each spanwise 
end of some of the chordwise modes. These node lines result from the 
rigid tip bulkheads used in the construction of the specimens. 

Since the covers of the monocoque beam were uniform in thickness 
and the camber of the cross section can be closely approximated by a 
parabola, the exact solution previously discussed is applicable. Hence 
the predominantly chordwise and predominantly spanwise frequencies of 
the monocoque including the influence of coupling can be obtained 
directly from figures 5 and 8. 
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Before this can be done, however, the uncoupled beam frequencies 
“r, B and uncoupled cross-sectional frequency a) c must 

be obtained. The uncoupled chordwise frequency *« was calculated 

-L 

with the use of the curves in figure 4(a) or equation (39) for symmetrical 
frequencies and equation (34) (with k = 0) for the antisymmetrical case. 
The uncoupled chordwise frequencies are listed at the top of tables I 
and II. The bending frequencies of the beam were calculated in refer- 
ence 2 by methods which include effects of transverse shear and shear 
lag, and the torsional beam frequencies by methods which include restraint 
of warping. These beam frequencies are listed in column (2) of tables I 
and II. 


In the first column of tables I and II is shown a sketch of the 
nodal pattern under consideration. The modes are grouped in pairs, 
except for the first, with each pair consisting of a spanwise mode and 
the corresponding first coupled chordwise mode. (The higher chordwise 
modes were not obtained experimentally and thus are not included. ) The 
first mode in column (T) of each table is the chordwise mode associated 


CO -| (JD T 

with each of the rigid- body spanwise modes. The ratios c 2 — and c * 1 




ox 


O) 


r, T 


are shown in column ( 
sponding frequency ratios 


of tables I and II, 
o> ^ 0) 


and 


r, B 

respectively. The corre- 
were obtained from figures 5 


“^B ' °V,T 

and 8 and are tabulated in column (4) . For the predominantly beam modes 
this ratio represents the reduction in beam frequency due to cross- 
sectional deformations. In an attempt to evaluate the results, the cal- 
culated coupled frequency co Is presented in column (3) and the experi- 
mental frequency in column © . 


In table I (for bending) the predominantly beam modes show little 
effect due to coupling at the low frequencies. At the higher modes, 
however, the coupling effect becomes appreciable, and it is encouraging 
to note the good agreement between calculated and experimental frequencies. 
This good agreement at the higher modes can also be looked upon as an 
indication that the method is applicable to lov-aspect-ratio wings. 

In table II (for torsion) the predominantly beam modes are greatly 
affected by the coupling due to cross-sectional deformations, even at 
the low frequencies. However, the accuracy of prediction of the coupled 
frequencies at the higher modes indicates that the applicability of the 
method in torsion to low- aspect- ratio wings has limitations. 


The calculated coupled chordwise frequencies show fair agreement 
with experimental frequencies and indicate that a large number of chord- 
vise modes are present over a narrow frequency band, sometimes falling 
within a few cycles of each other. Unfortunately, as the investigation 
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presented in reference 2 was primarily concerned with beam modes, the 
experimental frequencies of the chordwise modes were not always obtained 
and some are missing in tables I and II. Further, because of the com- 
plexity of the nodal patterns, it was difficult in some instances to 
correlate the experimental and theoretical mode shapes. Though the cor- 
relation between the calculated and experimental chordwise frequencies 
is good, the application to low-aspect-ratio wings is open to question, 
as the accuracy of prediction is not well established at the higher 
frequencies. 

It is to be noted that a second, thicker skinned monocoque beam was L 

discussed in reference 2, but for a number 

sented or discussed here. First, the f — 

\ho 

thicker monocoque is felt to be straining the limits of some of the basic 
assumptions of the present analysis. Second, because of a much higher 
uncoupled chordwise frequency, the coupling effects on the beam frequencies 
in bending would be much smaller than in the lighter monocoque. Further- 
more, the chord length and radius of curvature used in the calculations 
of the uncoupled beam frequencies g and o> r ^ are questionable, 

as the frequencies of the rigid- body beam modes g calculated from 

these dimensions are lower than those observed experimentally. Since 
the theories used in these calculations do not include all the flexibil- 
ities of the actual structure, the calculated values should be higher 
than the experimental values, rather than lower. It may also be noted 
that the dimensions of maximum depth, chord length, and radius of curva- 
ture for the thicker beam are not compatible with those for a circular 
arc. When the radius of curvature was altered by 6 percent in order to 
describe a circular arc, a 12-percent change was noted in the third and 
fifth rigid- body modes. More than likely, then, the calculated rigid- 
body frequencies of the thicker beam are somewhat higher in magnitude 
than those calculated in reference 2. 


of reasons it has not been pre- 
2 

ratio of 0.356 for the 
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CONCLUDING REMARKS 


The vibrational behavior of thin wings with airfoil-type cross 
sections has been investigated by means of an analysis of an idealized 
structure consisting of two curved beams mounted on springs. The results 
obtained from an analysis of this model can be* used not only to obtain 
the chordwise modes and frequencies of the prototype but also to obtain 
an estimate of the reduction in spanwise frequency due to the coupling 
effect between cross-sectional and spanwise motions. Both uncoupled 
chordwise frequencies and the resulting coupled chordwise and spanwise 
frequencies are presented for beams with uniform cover thickness and 
parabolic cross sections. A comparison is made with existing experimental 
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data on a circular-arc monocoque. The comparison shows good agreement 
for the coupled beam frequencies and fair agreement for the coupled 
chordwise frequencies at moderate aspect ratios. For beams of arbitrary 
cross section an approximate matrix procedure is also presented. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Air Force Base, Va. , September 21, 1961 . 
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APPENDIX 

APPROXIMATE SOLUTION OF AN ARBITRARY AIRFOIL SECTION 


In this appendix an approximate solution for an idealized structure 
of arbitrary shape and variable thickness is presented. The energies of 
the system are expressed in discrete form by writing all derivatives as 
differences and all integrals as summations according to the trapezoidal 
rule. The deflections at discrete points are considered as the unknowns L 

of the system. After the boundary conditions are properly accounted for, 1 

a system of linear simultaneous equations in terms of these unknowns is 4 

obtained by minimizing the energy function U - T. ^ 

4 

From equations ( 4 ), (6), and (7) the expression for U - T can be 
written as 

U - T = - 
2 




The first step in the procedure is to choose equally spaced stations 
along the x-axis, numbering them from n = 0 to n = N. Integrals can 
then be approximated, either in terms of values of the integrand at the 
full stations 


f(x)dx = + f x + . • . + ^ (A2) 

or in terms of the values at the half stations 


f 0 f(x,dX ' A [ f l/2 + f 3 /2 + ' ' • + V(l/2)] <«> 

where A is the equal distance between stations. 

The first integral in equation (Al) is approximated by means of 
equation (A2), as the values of the second derivative appearing in the 
integrand are best approximated by differences at the full stations: 



25 


L 

1 

4 

4 

4 


w 


w - 2w + w , 
n+1 n n-1 


(A4) 


Consequently, the first integral can be written in matrix form as follows 


El(w") 2 dx = ^L W J[V]'[ EI ][^*]| W 


(A5) 


where the definitions of the matrices are given below, and the number of 
rows and columns, in that order, are indicated beside each: 



W N+1 


(!)x(N+3) (a6) 


W 



(N+5)x(l) (A7) 



(n+i)x(n+3) (a8) 


1-2 1 
1-2 1 



(n+i)x(n+i) (A9) 


A prime on a matrix denotes the transpose of the matrix. 
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Note that the values of w at stations corresponding to n = -1 
and n = N + 1 (points not on the structure) appear in the matrices. 
These points can be evaluated with the use of the finite-difference form 
of the boundary conditions in equations (10) and (12): 


and 


v. 


EIv" 


n+1 


- 2w n + v- 


N-l 


I n=N 


= 0 


EIw" 


Vi - 2 w q + w_ x 


n=0 


= 0 


(A10) 


(All) 


With the conditions expressed in equations (A10) and (All), equa- 
tion (A5) can be rewritten as 


L 

1 

4 

4 

4 


I f ° Ei(u'') 2 dx - ^[*J[d 2 ] , [ei][d 2 ]|w| (A12) 


where now 




1-2 1 
1-2 1 

1-2 1 
1-2 1 


(1)X(N+1) (A13) 


(n-i)x(n+i) (Al4) 


r. 


M ■ 


(EI). 


(El) 


N-] 


(N-l)x(N-l) (A15) 
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The second integral in equation (Al) is approximated by the use of 
the values at the half stations (eq. (A3)), as the value of the first 
derivative appearing in the integrand is, in turn, best approximated by 
a half-station expression: 


v 


_ v n+l “ w n 

n+(l/2) * A 


(Al 6) 


Thus, by using both the half-station expressions (eqs. (A3) and (Al6)) 
the second integral can be written in terns of the displacements at the 
full stations : 


i^EthVaxjf h'wdx -|L w J[ D l] , [h'jM[h(|[l^|v| (A17) 


Here 



(l)x(N+l) 


(A18) 



(N+1)X(1) (A19) 


-1 


1 




h 3/2 


-1 


1 


-1 1 


-1 


(n)x(n+i) 


(A20) 


h N-(3/2) 

h N-(l/2) 


(N)x(N) (A21) 



2 6 


(Et) l/2 

(Et) l/2 

(Et)j/ 2 

<^>5/2 


^ Et ^N— (l/2) ^^-(1/2) 


(N)x(N) (A22) 


The third term in equation (Al), or the potential energy of the 
spring mounts^ is given in matrix form by 


i k^c) 2 + i k 2 w(0) 2 



(A23) 
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where 




(N+l)x(N+l) 


(A24) 


The last integral in equation (Al) is approximated by equation (A2): 



(A25) 


where 




(N+l)x(N+l) (A26) 
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(eq. (S) e Ste®!“l’ M to). (A2J), and (A2 5 ) 
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the 


where 


2 (Bjliv 


energy function 

(A27) 


W • IH'hh MW . jj,j m 

[ B j = AjTptj 


S3“ n=rs: ;e~ “ r:t s . 


(A29) 


Equation (A30) can be 


since 


jf A M A ]'] -T!W + [Bj']j) w j = 0 

rewritten in a more usable form as 

w| ft'fw - wjv| 

m* ■ w 


(A30) 


(A31) 


(A32) 


and the transpose of a i s 
1 ~^ ord. 

toaohdeva 

calculated from the e S bectlon shape of the first svL a b f ief comparison 
symmetrical parabolic Ct solutlon Presented for the^T" 6 * 1 ' 108 ' 1 m ° de as 

comparison, only half ZTIu S ^ nCe a symmetrical mode Vfi : J * d by 

1 ’ hUS beC °" e 0 “1 c /2 instead^of 
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. , _, 1fl cation of constant cover thick- 

0-d c.) equation (A 5 l) reduces to 

ness and constant materiai P v 

criptu) 


with 


[b] “ 


(N+l)x(N+l) 


(A55) 


M'S? 


(N+1)X(N+1) (A56) 


(A3T) 


A = i n 


-.«• .. - J-M S WfJS”.*? 

sssmtsst* 1 


W -1 = W 1 


end Dg in equation (A53) ^comes 


D 3 


ir. equation ( A3& ) • 



1 


-2 1 

1 -2 1 I 
1-2 1 


(N)X(N+1) (A39) 
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P' -P' 4=~ h tr 1 
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Furthermore, 



1 

1 


(N)x(N) 


v* 

Jl 

2 




5 



(N)x(N) 


L 



T) 


-^1 



(N)x(N) 




1 

-1 1 


-1 1 
-1 1 
-1 1 


(n)x(n+i) 


0 



0 


L 


0 


il 


(N+l)x(N+l) 


(A40) 


(Akl) 


(Ak2) 


(Akl) 


(Akk) 
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If the rigid-body (translational) frequency ie represented as 


2 = JL. 

“r, B cpt 


with 



4 _ c^pta)^ 

a “ i6ei 


equation (A34) can be expressed in the parameters of the exact solution 
and becomes 



In figure 10 a plot of the first symmetrical mode shape of the 

ir - n nn(i (— \ = 0.018 is shown. The 

parabolic -arc monocoque for k - O ana f j 
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Figure 5>- Concluded. 
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